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We report a measurement of the branching fraction of the decay B° — > D* + D*~ and of the CP-odd 
component of its final state using the BABAR detector. With data corresponding to an integrated 
luminosity of 20.4 fb" 1 collected at the T(4S) resonance during 1999-2000, we have reconstructed 38 
candidate signal events in the mode B° —> D* + D*~ with an estimated background of 6.2±0.5 events. 
From these events, we determine the branching fraction to be B(B° —> D* + D*~) — (8.3± 1.6(stat) ± 
1.2(syst)) x 10" 4 . The measured CP-odd fraction of the final state is 0.22 ± 0A8{stat) ± 0.03{syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



After the observation of time-dependent CP-violating 
asymmetries in the decays of neutral B mesons to CP 
eigenstates containing charmonium Q it is inter- 
esting to extend the search for CP-violating effects to 



Cabibbo-suppressed double charm modes, such as B — > 
£)(*)+£)(*)- ||. The interference of the dominant tree 
amplitude with the mixing diagram is sensitive to the 
angle [3 of the Unitarity Triangle in this case as well; 
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however, the theoretically uncertain contribution of pen- 
guin amplitudes with different weak phases is potentially 
significant and may shift the observed asymmetry by an 
amount that depends on the ratio of the penguin and 
tree contributions and their relative weak phases. The 
B° — > D* + D*~ vector-vector final state has very clear 
experimental signatures that make it an interesting can- 
didate for CP-violation measurements. However, it is 
not a pure CP eigenstate and a dilution of the measured 
asymmetry can be produced by a P-wave, CP-odd, com- 
ponent. A time-dependent angular analysis of the de- 
cay products pj can remove the dilution by resolving 
the CP-even and CP-odd components. As a precursor 
to measuring time-dependent CP-violating asymmetries 
using the decay B° — > D* + D*~ , we report in this letter 
a measurement of the B° — > D* + D*~ branching frac- 
tion and a measurement of the CP-odd component, R±, 
of the final state. These measurements represent sig- 
nificant improvements over the previous measurements 
B(B° -> D*+D*-) = (9.9tH(stat) ± 1.2(syst)) x 1CT 4 
and (1 - R ± ) < 0.11 at 90%'C.L. [§. 

The data used in this analysis were collected with the 
BABAR detector [g at the PEP-II storage ring @ located 
at the Stanford Linear Accelerator Center. This data 
sample represents an integrated luminosity of 20.4 fb _1 
collected on the T(AS) resonance. Assuming 50% of the 
T(4S) decays give B°B°, the number of neutral B mesons 
in this sample is (22.8 ± 0.4) x 10 6 . 

Charged particles are detected and their momenta 
measured with the combination of a 40-layer drift cham- 
ber (DCH) and a five-layer silicon vertex tracker (SVT), 
both operating in a 1.5 T solenoidal magnetic field. The 
charged particle tracking system allows particles with low 
momentum in the laboratory frame to be reconstructed 
efficiently, a property that is very important for this anal- 
ysis. This efficiency begins to turn on at a momentum 
of ~60 MeV/c and reaches its maximum value at around 
200MeV/c. Photons are detected by a CsI(Tl) electro- 
magnetic calorimeter (EMC) that provides high detec- 
tion for energies above 20MeV, with typical energy and 
angular resolutions of 3% and 4mrad, respectively, for 
1 GeV photons. Charged particle identification is pro- 
vided by the ionization loss measurements in the SVT 
and DCH, and by an internally reflecting ring-imaging 
Cherenkov detector (DIRC) covering the central region 
of the detector. 

Events are selected by requiring three or more charged 
tracks and that the normalized second Fox- Wolfram mo- 
ment H of the event be less than 0.6. We also re- 
quire that the cosine of the angle between the recon- 
structed B direction and the thrust axis of the rest of 
the event, calculated in the T(AS) rest frame, be less 
than 0.9. These criteria discriminate T{AS) events from 
non-resonant background events. 

B° mesons are exclusively reconstructed by combin- 
ing two charged D* candidates, using a number of 



D* + and D decay modes. The D* + mesons are re- 
constructed in their decays D* + — > D°n + and D* + —> 
D + tt°. We include in this analysis the decay combi- 
nations P>* + P>*-decaying to (D°tt+, D° tt") or (P>°7r+, 
D~tt°), but not (D + tt° ,D~tt°) due to the smaller branch- 
ing fraction and larger expected backgrounds. D and 
D + candidates are subjected to a mass-constrained fit 
to provide an improved measurement of the D meson's 
momentum. They are combined with pion candidates, re- 
ferred to as "soft" pions due to their low (< 200MeV/c) 
transverse momentum, to form D* + candidates. A topo- 
logical vertex fit is performed that includes the mean 
position of the e + e~ interaction point to improve the an- 
gular resolution of the soft pion. 

The decay modes of the D° and D + are selected by an 
optimization of S 2 / (S + B) based on Monte Carlo simu- 
lations, where S and B are the expected number of signal 
and background events, repsectively. We first determine, 
based on Monte Carlo simulations, the expected S — B 
for each of the decay mode combinations individually. 
Then, we successively add modes in order of decreas- 
ing S — B to compute an overall S 2 / (S + B) value until 
S 2 /(S + B) no longer increases. The decay modes used 
are D° -> K-*k + , D° -> R-ir+ir , D° -> K-ir+Tr-ir+, 
D° -► K° s ir-ir+, D+ -> K--K+TI+, D+ -> K°n+, and 
D + — ► K~ K + tt + . D° (D + ) meson candidates are re- 
quired to have a reconstructed invariant mass within 
20MeV/c 2 of the nominal D° (£>+) mass {§. 

Charged kaon candidates are required to be incon- 
sistent with the pion hypothesis, as inferred from the 
Cherenkov angle measured by the DIRC and the ion- 
izations measured by the SVT and DCH. No particle 
identification is required for the kaon from the decay 
D° -> K--K+. 

K® — > 7r + 7r~ candidates are required to have an in- 
variant mass within 15MeV/c 2 of the nominal K® mass. 
The angle between the flight direction and the momen- 
tum vector of the K® candidate is required to be less 
than 200mrad, and the transverse flight distance from 
the primary event vertex, obtained from the remaining 
charged tracks in the event, must be greater than 2 mm. 

Neutral pion candidates are formed from pairs of pho- 
tons in the EMC with energy above 30MeV, an in- 
variant mass within 35MeV/c 2 of the nominal tt° mass, 
and a summed energy greater than 200 MeV. A mass- 
constrained fit is applied to these tt° candidates. In the 
case of the soft ir° from D* + — > D + ir° decays, the energy 
cut is replaced by a momentum cut, in the T(45 l ) frame, 
of 70 < p* < 450 MeV/c. 

To select B° candidates with well reconstructed D* + 
and D mesons, we form a x 2 that includes all measured 
D* + and D masses: 
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Anifl. — Amp* 



HID 
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Amjj* 



(1) 



where the caret over a value refers to the nominal value, 
and Amu* is the D* + — D mass difference. For a mn 
we use values computed for each D candidate, while for 
CAm D . we use fixed values of 0.83MeV/c 2 for D* + — > 
L»°7r+ and 1.18MeV/c 2 for £>*+ -> D+n". A requirement 
that xifass < 20 ^ s a PP ne d to all B° candidates. In 
events with more than one B° candidate, we choose the 
candidate with the lowest value of Xmoss ■ 

A B meson candidate is characterized by two kinematic 
variables: the energy-substituted mass, 



m ES = \/£f eam 



Pb 



(2) 



and the difference of the B candidate's energy from the 
beam energy, 



AE = E 



B 



E 



Beam ' 



(3) 



Eg {p*b) are the energy (momentum) of the B candi- 
date in the center-of-mass frame and E* Beam is one-half 
of the total center-of-mass energy. The signal region in 
the AE vs. m ES plane is defined to be \AE\ < 25MeV 
and 5.273 < to es < 5.285 GeV/c 2 . Based on Monte Carlo 
simulations, the width of this region corresponds to ap- 
proximately ±2.5cr in both AE and Toes- 

To determine the expected contribution from back- 
ground in the signal region, we scale the number of events 
seen in a sideband in the AE vs. jjjes plane defined as 
| AS | < 200 MeV, 5.20GeV/c 2 < to es < 5.26GeV/c 2 and 
50MeV < \AE\ < 200 MeV, 5.26 < m ES < 5.29GeV/c 2 . 
The scaling factor is calculated by parameterizing the 
shape of the background in the AE vs. m E s plane as the 
product of an ARGUS function in m E s and a first- 
order polynomial in AE . Based on this parameterization 
we estimate that the ratio of the number of background 
events in the signal region to the number of events in the 
sideband region is (1.72 ± 0.10) x 10~ 2 . The uncertainty 
is derived from the observed variation of this ratio un- 
der alternative assumptions for the background shape in 
m E s and AE using Monte Carlo simulations. The simu- 
lations also indicate that there are no significant sources 
of background appearing in the signal region beyond that 
indicated by the sideband extrapolation. 

After all selection criteria, 38 events are located in the 
signal region, with 363 events in the sideband region. The 
latter, together with the scaling factor determined above, 
implies an expected number of background events in the 
signal region of 6.24±0.33(stai)±0.36(sys£). The system- 
atic uncertainty comes from the background shape varia- 
tion mentioned previously. Figure [l] shows a projection of 
the data on the m E s axis after requiring \AE\ < 25 MeV. 
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FIG. 1: The m E s istribution of 5° -> D* + D*~~ events with 
\AE\ < 25 MeV. The curve represents a fit with the sum of 
a Gaussian to model the signal and an ARGUS function jnj] 
to model the background shape. 



We use a Monte Carlo simulation of the BABAR de- 
tector to determine the efficiency for reconstructing the 
signal. The efficiencies range from 17.4% to 2.7%, de- 
pending on the D decay modes. This, together with the 
total number of neutral B mesons produced during data 
collection, allows us to determine the branching fraction 
for B° -> D*+D*~ to be 



B(B° 



D*+D* 



(8.3±1.6(stat)±1.2{syst)) x 10" 



The high charged particle multiplicity makes this mea- 
surement particularly sensitive to the tracking system. 
Therefore the dominant systematic uncertainty comes 
from our level of understanding of the charged particle 
tracking efficiency. Systematic errors are assigned on a 
per track basis for n, K, and soft it, and are added lin- 
early (9.9%). The effect on acceptance due to the impre- 
cisely known partial-wave content of the B a — > D* + D*~ 
final state is another source of potential systematic bias 
(6.6%). Other significant potential systematic biases 
arise due to the uncertainties on the branching frac- 
tions § of the £>* + , D°, and D+ (5.6%) and the un- 
certainties in mass resolutions of reconstructed mesons 
(4.1%). The total systematic uncertainty from all con- 
sidered sources is 15%. 

In addition to the branching fraction quoted above, 
we have also measured the CP-odd fraction of the final 
state. This fraction, R± , is determined from the angular 
distribution of the soft pions in the decay, analyzed in 
the transversity basis B. In this reference frame, three 
decay amplitudes determine the distribution of three de- 
cay angles. Integrating over time, B flavor, and two of 
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these three angles yields the following expression: 



1 oT 



T d COS t r 



-(l-R±) sin 2 9 tr + t;R± cos 2 



(4) 



Here T is the decay rate and 9 tr is the angle between the 
normal to the D* ~ decay plane and the line of flight of 
the soft pion from the D* + evaluated in the D* + rest 
frame. 

We perform an unbinned maximum likelihood (ML) 
fit to the 38 events in the signal region described previ- 
ously. The fit takes into account the presence of back- 
ground, whose properties are derived from the sideband 
sample, and the angular resolution ag estimated from 
Monte Carlo simulations. We define the likelihood func- 
tion to be 



i— l,n 

(l-p)xF(6tr t i,a8,i,R"x g ) > 



i—l,n i—l.n 

bkg\ 



where n is the number of selected events. The contri- 
bution to the total likelihood from the z-th event, Ci 1 
is defined in terms of the purity, p, of the sample and 
the probability density functions J-(9t r ,i, &e,i, R±) for the 
signal and background. R^ 9 and R± 3 are the parame- 
ters describing the shapes of the signal and background 
angular distributions, respectively, and 6t r ,i is the mea- 
sured transversity angle in event i. The probability den- 
sity functions T are obtained from the convolution of 
the angular distribution (Eq. ^J) with Gaussian resolu- 
tion functions describing the measurement uncertainties 
ag t i. From studies of simulated data, ag was measured 
to be 0.11 (0.12) radians for charged (neutral) slow pi- 
ons. A 10% uncertainty on these values is considered 
when estimating the corresponding systematic error. 

The value of R b ± 9 is evaluated by fitting the 363 events 
in the sideband region and setting p = in Eq. |^. The 
result of this fit is R b ± 9 = 0.29±0.04, compatible with the 
value expected for a uniform distribution (R± = 1/3). To 
determine R s , g , we fit the 38 events in the signal region 
with i?f 9 fixed to 0.29 and p fixed at 83.6%. The result 
of the fit to the signal region, without the correction for 
angular acceptance bias described below, is R 8 ^ 9 — 0.25 ± 
0.18(stat), and is shown in Fig. |[ The probability of 
obtaining a lower likelihood, evaluated using a Monte 
Carlo technique, is 66%. 

It should be noted that Eq. [| is the differential decay 
rate T integrated over the full ranges of the other two de- 
cay angles in the transversity basis, neglecting any bias in 
the projected tr distribution introduced by detector ac- 
ceptance effects. A detailed study of the kinematics of the 
decay shows that the incomplete detector coverage of the 
polar angle with respect to the beam axis does not intro- 
duce any bias in the distributions of the decay angles in 
the transversity basis. However, an inefficiency in detect- 
ing soft pions below a threshold in transverse momentum 
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FIG. 2: The cos6* tr distribution from the unbinned ML fit, 
superimposed on the histogram of the B° — > D* + D*~ candi- 
dates in the signal region. The solid line represents the cos 6 tr 
distribution from the unbinned ML fit for the selected events. 
The dotted and dashed lines represent the fitted CP-odd and 
CP-even components, respectively, for the signal. Lhe dot- 
dashed curve represents the fitted background component. 



may indeed introduce such a bias due to the correlations 
between the decay angles and particle momenta in the 
laboratory frame. The amount of these acceptance losses 
depends on the population of phase space, determined by 
the values of the decay amplitudes. 

An accurate correction for these acceptance effects re- 
quires the complete determination of the decay ampli- 
tudes using a full angular analysis on a sufficiently large 
data sample. To estimate the size of the acceptance bias 
on R± without knowing the decay amplitudes, the fit pro- 
cedure was tested on several samples oi B° — > D* + D*~ 
simulated events generated with different sets of decay 
amplitudes. The different amplitudes affect to varying 
extents the correlated soft pions' transverse momenta and 
angular distributions. The fitted R± values were found 
to be consistent with the generated values in the limit of 
negligible soft pion inefficiency. Depending on the mix of 
decay amplitudes, they did reveal a bias once the pion- 
detcction threshold was taken into account. Considering 
the full possible range of decay amplitudes, the calculated 
bias on R± ranged from -0.048 to +0.004. The central 
value of this interval is taken as a correction to the fit- 
ted i?^ 9 , while its half width is taken as an estimate of 
the corresponding systematic uncertainty (0.026). Addi- 
tional, smaller systematic uncertainties affecting the R± 
measurement arise from the imperfect knowledge of the 
resolution in the transversity angle &tr (0.006), the angu- 
lar distribution of the background (0.008), and the purity 
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of the signal sample (0.0003). The total systematic un- 
certainty on R± is determined to be 0.028, giving the 
final corrected result: 

R ± = 0.22 ± 0.18(stat) ± 0.03(syst). 

In summary, we have observed a signal of 31.8 ± 
6.2(stat) ±0A(syst) events in the decay B° -> D*+D*-. 
Our measurement of the branching ratio is 

B(B° -> D*+D*-) = {8.3±l.6{stat)±1.2(syst)) x 10~ 4 . 

From the transversity angular distribution of these 
events, we have also measured the CP-odd fraction, R±, 
of the final state. These measurements provide a starting 
point for measuring time-dependent CP-violating asym- 
metries in these decays when more data become available. 
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